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Localized single-particle potentials for all octet baryons, A^, A, S, and H, in finite nuclei, '^^C, 
16q^ 28gj^ *°Ca, ^^Fe, and ^°Zr, are calculated using the quark-model baryon-baryon interactions. G- 
matrices evaluated in symmetric nuclear matter in the lowest order Brueckner theory are applied to 
finite nuclei in local density approximation. Non-local potentials are localized by a zero-momentum 
Wigner transformation. Empirical single-particle properties of the nucleon and the A hyperon in 
nuclear medium have been known to be explained semi-quantitatively in the LOBT framework. 
Attention is focused in the present consideration on predictions for the E and H hyperons. The 
unified description for the octet baryon-baryon interactions by the SUe quark-model enables us to 
obtain less ambiguous extrapolation to the S = —1 and S = —2 sectors based on the knowledge 
in the A'^A'^ sector than other potential models. The E mean field is shown to be weakly attractive 
at the surface, but turns to be repulsive inside, which is consistent with the experimental evidence. 
The H hyperon s.p. potential is also attractive at the nuclear surface region, and inside fiuctuates 
around zero. Hence H hypernuclear bound states are unlikely. We also evaluate energy shifts of the 
E~ and H~ atomic levels in ^*Si and ^'^Fe, using the calculated s.p. potentials. 



PACS numbers: 21.30.Fe, 21.65.-f, 21.80.+a 

I. INTRODUCTION 

One of the salient features of atomic nuclei is the suc- 
cess of the description of their properties by a single- 
particle (s.p.) picture. Although the nucleon-nucleon 
interaction is known to be strongly repulsive at the 
short range part, which was once conveniently described 
by a hard-core, the nucleon single-particle potential is 
well represented by a well-behaved local potential of the 
Woods-Saxon form. The theoretical base of understand- 
ing this circumstance in view of the singular two-body in- 
teraction was provided by the Brueckner theory in 1950's 
[U, 0, 3 . The progress of the density-dependent Hartree- 
Fock (DDHF) description of nuclear bulk properties fol- 
lowed in 1970's 0,01, introducing some phenomenological 
adjustment for G-matrices in the Brueckner theory. 

The mean field picture seems to hold also for hyper- 
ons in nuclei. For the A hyperon, the potential proper- 
ties have been known from light to heavy nuclei from A 
formation and spectroscopy experiments [6| . Experimen- 
tal studies of the E and S hyperons in nuclear medium 
properties are now in progress. Because direct hyperon- 
nucleon scattering experiments are not readily available, 
the properties of the hyperon s.p. potentials are a valu- 
able source of hyperon-nucleon interactions. This case, 
we have to resort to an effective interaction theory to re- 
late s.p. properties of the hyperon embedded in nuclei 
with the character of hyperon-nucleon interactions. 

In this paper we develop a method to obtain local po- 
tentials for octet baryons in finite nuclei with using full 
non-local G-matrix elements in nuclear matter, starting 
from the baryon-baryon bare interactions. The calcula- 



tion of single-particle properties in nuclear matter can 
provide the basic information about the baryons in nu- 
clear medium derived from the bare interaction. Nev- 
ertheless, it is instructive to explicitly calculate the s.p. 
potential in finite nuclei starting from two-body baryon- 
baryon interactions and compare them with the empirical 
ones. It is not even obvious whether the shape repre- 
sented by the Woods-Saxon form which has been estab- 
lished both for the nucleon and the A hyperon mean fields 
is suitable for the E and the S hyperons. The straightfor- 
ward folding of the two-body effective interaction in mo- 
mentum space provides a non-local potential in a nucleus. 
The non-locality also comes from the exchange character 
of the basic interaction. In order to make a comparison 
with empirical data, it is meaningful to define a local po- 
tential by some localization procedure. We employ in this 
paper a zero-momentum Wigner transformation method 
based on the WKB localization approximation 7]. 

It is necessary for a predictive discussion about hy- 
peron s.p. properties in nuclei to use octet baryon-baryon 
bare interactions as reliable as possible. With little ex- 
perimental information except for the AA'^ interaction, 
the construction of the interactions in the strangeness 
S — —1 and S = —2 sectors is not simple, although 
some constraints are imposed by the flavor symmetry. 
The typical potential model has been developed in a 
one-boson-exchange potential (OBEP) picture by the Ni- 
jmegen group. The early parametrization with the hard- 
core in the 1970s [sl, 0] has been successively revised by 
adjusting parameters in the soft core version pH, [ll| and 
introducing new terms p^ . There are now a number of 
sets of parameters, reflecting ambiguities due to the lack 
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of experimental data. Although the description for the 
AN interaction seems to be under control, there are var- 
ious uncertainties in the 'SN sector. For the S hyperon, 
namely the strangeness S* = — 2 sector, the situation is 
no better. 

Using the s pin- flavor SUe quark-model, the Kyoto- 
Niigata group m, [H, m m have developed a uni- 
fied description for the interactions between full octet 
baryons. In this model, the interaction is constructed 
as the Born kernel in the framework of the resonating 
group method for the three-quark clusters, the short 
range part of which is composed of the effective one-gluon 
exchange mechanism. The basic SUg symmetry provides 
a specific framework to the interactions between octet 
baryons and the Pauli principle respected on the quark 
level in addition brings about a characteristic structure to 
them. Incorporating effective meson exchange potentials 
between quarks, namely the scalar, pseudo-scalar, and 
vector mesons exchanges, the model is able to account 
the NN scattering data as accurately as other modern 
NN potential models. 

Parameters of the SUg quark-model by the Kyoto- 
Niigata group are mostly fixed in the NN and AA^ sec- 
tors, and the uncertainties in the extension to the "SN 
and EN channels are limited. In fact, the definite pre- 
dictions such as the smallness of the AA^ spin-orbit inter- 
action and the overall repulsive nature of the E-nucleus 
s.p. potential in nuclei have been supported by the exper- 
iments afterwords. Therefore it is interesting to examine 
the whole prediction of this potential for s.p. properties 
of all the octet baryons in nuclear medium. In particular, 
concrete predictions are presented for the S hyperon. We 
use the most recent quark-model potential fss2 
in this paper. The parameter set includes no adjustable 
parameter for the tuning afterward. The original interac- 
tion as a Born kernel has an inherent energy-dependence. 
Recently the method to eliminate the energy-dependence 
was developed [3]. We actually use this renormalized 
version of the fss2 potential. 

We present, in Sec. II, basic expressions of the method 
to evaluate localized N, A, E, and S s.p. potentials in 
a finite nucleus. We first discuss numerical results of 
these s.p. potentials in nuclear matter to represent basic 
characters of the G-matrices of the quark-model baryon- 
baryon interactions. Calculated results in finite nuclei, 
12c, 1^0, 28si, 40qj^^ 56pg^ 90^^^ ^rc shown in Sec. 
IV. The energy shift and the width of the and 



atomic levels in ^^Si and ^^Fe are studied in Sec. V on the 
basis of the s.p. potential obtained in Sec. IV. Summary 
is given in Sec. VI. 



II. LOCALIZED SINGLE-PARTICLE 
POTENTIAL IN A FINITE NUCLEUS 

We calculate a baryon single-particle potential in a fi- 
nite nucleus which is defined by folding A^A^ or YN G- 
matrix elements in nuclear matter with respect to nu- 
cleon occupied states through the local density approx- 
imation. It has been a traditional method for the mi- 
croscopic study of bulk properties of nuclei to construct 
density-dependent two-body local interaction based on 
the G-matrices in nuclear matter and apply the effec- 
tive interaction to mean field calculations for finite nu- 
clei. Avoiding this procedure, we directly fold G-matrix 
elements to obtain non-local s.p. potentials and localize 
them. In the DDHF calculations, some phenomenological 
adjustments are introduced to reproduce the properties 
of the well-known nuclei. The purpose of the present 
paper is not to accomplish the reproduction of empiri- 
cal values, but examine overall implications of the uni- 
fied description of the baryon-baryon bare interaction 
by the quark model potential fss2 for hyperon s.p. 
potentials in finite nuclei. In this section, we derive a 
basic expression for the s.p. potential, by introducing 
some approximations and a localized method by the zero- 
momentum Wigner transform. 



A. direct term 

First we consider the following direct term contribu- 
tion. The wave function (jjif^j^mj^ denotes the nucleon 
s.p. wave function of the nucleus with the orbital angu- 
lar momentum £h and the total angular momentum jh, 
and the dummy wave function of the baryon for which 
the potential is calculated is denoted by 4'ijmj ■ The av- 
erage over the z-component of the total angular momen- 
tum, J with j = 2j + 1, means that we assume 
the spherical symmetry from the beginning. We do not 
write the isospin indices for simplicity in the following 
derivation, but recover them in the final expression. 



^ JJJJ dridr2dr\dr'2(j}*gj„^^{r[)(j)}^j^^,^^^{r'2)G{r\,r'.2,ri,r2)^^^ 

r £ £h L] ( £ 4 ^' 1 . . . . 
= ^ ^ jhSVLU I 1/2 1/2 5 W 1/2 1/2 S \ dridr2dr[dr'2 

h JMLL'S [ j jh J ) [ j jh J ]•'•'•' •' 
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The effective baryon-baryon interaction G{r[,r'2,ri,r2) in a coordinate space is supposed to be related to the G- 
matrix in momentum space G{k\k; K,lu) by 

GK,r^,ri,r2) = JJJJ dk,dk,dk[dKSiK - K')e^^'''^-'''^+k^^^^^^ (2) 

Here each momentum has the following relation: fci = — — K + k, k2 = — — K — k, k\ = — 2|j — x' + fc', 

^ ^ mi+m2 ' m\+m2 ' -I- mi+m2 ' 

k'2 — „-^"^^-^^ K' — k' . The mass of the baryon for which the s.p. potential is evaluated is denoted by toi and the mass 
of the nucleon in the target nucleus by 1112 ■ The G-matrix is evaluated in symmetric nuclear matter by solving the 
baryon-channel coupling Bethe-Goldstone equations 

G^,^,{k',k;K,u;)^V^,^,{k\k;K) + -^y2 I dqVo^Ak' K) ^^^'^i , G/3,a' (g, fc; i^, (3) 

K^'^) a J - Eb{ki) - EN{k2) 



with the suffix specifying the pair of a baryon b and a nu- 
cleon N by p. The Pauli exclusion operator Qfj is treated 
in the standard angle-average approximation. The ex- 
plicit expression may be found in ref. (20| . Ea{k) is a 
s.p. energy of the baryon a in nuclear matter. We em- 
ploy the continuous choice for the energy denominator 
in Eq. (3). That is, Ea{k) ^ rua + + Ua{k) is 

defined self-consistently by the following definition of the 
s.p. potential Ua- 



Ua{k) = J dk'Ga^a{q,q;k + k',uj), 



(4) 



where q = \{k-k') and uj = EN,Y{k) + E^ik'). The 
prescription for the starting energy lj in the local density 
approximation is explained in the next section. 

The straightforward calculation of Eq. (1) needs much 
computational effort and is not fruitful to obtain a phys- 
ical insight for baryon properties in nuclei by starting 
from the bare baryon-baryon interactions. We introduce 
two simplifying approximations. One is the spin-average 
in taking the sum of the matrix elements, which means 
that we take the average over the spin orientation: 

iY.^SMs\G{k',k-K,u)\SMs) 

^ Ms 

qJq ^ 

where the G-matrix is decomposed to partial waves and 
Pq stands for the Legendre polynomial with q specifying 
the orbital angular momentum. The other simplification 
is the following replacement. 



dRidsict)*^^^. [Ri + ]-si)(j)ij^. {Ri ~ ^si) 



=^ j dRi(j)*fj^^{Ri)4>ijm^[Ri) j dsi. 



(6) 



This corresponds to the zero-momentum Wigner trans- 
formation of the non-local potential. That is, we set 



p = for the Wigner transformation U^{R,p) of the 
non-local potential [/(ri,r2). 



:R,p)^ J dse'P-^U{R+^s,R-^. 



(7) 



Results shown in Sect. IV for nucleons and lambdas, for 
which we know what s.p. potentials potentials are ex- 
pected in G-matrix calculations with bare NN and A7V 
interactions in various studies in literature, implies that 
the zero-momentum approximation works well. More di- 
rect confirmation of the reliability of this approximation 
will be presented elsewhere (2l|. 

To evaluate Eq. (1) with the above simplification it is 
convenient to use the Fourier transform of the s.p. wave 
function 



(fc) 



1 



(27r)3 
1 
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dre-''^-^(t>ir) 



1 ~ 



i'--'[Ye{k)xxi/2Ym,i^iiik)l8) 



where n is a nodal quantum number and the Fourier 
transformation of the radial wave function is defined as 



dr rji{kr)(j)ij{r) 



(9) 



After carrying out some integrations and taking angular 
momentum recouplings, we obtain the final expression as 
follows. 



In = 



X / / dkdk'jo{\k' - k\Ri) 
1 



X 



IQil 



Q'll 

,(|Qi|)P,,(cosQ^) 



.JIQ'il) 



X E -^^G^'f (fc,A:')^,(cosfcfc'), 

qJ,S 



(10) 
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where Qi and Q[ are defined by k and fe' as 

B. Exchange term 

We also have to consider the exchange term contribu- 
tion, which is famihar for the nucleon through the an- 
tisymmctrization of the wave function. For hyperons, 
such terms appear in association with the exchange char- 
acter of the hyperon-nucleon interaction, which is real- 
ized by the strange meson exchange in the OBEP de- 
scription. Denoting the space-exchange operator and the 
spin-exchange operator by Pr and P^, respectively, and 
specifying the even and odd components of the interac- 



tion under the space-exchange, the matrix element of the 
YN interaction is written as 

{YN\V\YN) = {YN\VE^{l + Pr) + Vol{l-Pr)\YN) 
= {YN\^{Ve + Vo)\YN) 

-{YN\^{Vo-Ve)P.\NY), (13) 

where the relation P^Pr\NY) = \YN) is used. The first 
term was treated in the previous subsection as a direct 
term contribution, and the second term is considered in 
this subsection. The effective interactions in the direct 
and exchange contributions should be treated as such 
a combination of the even and odd parts in each spin 
and isospin channels, though the isospin dependence is 
disregarded in the above expression because the inclusion 
of it in the final expression is simple. 



rrij hrrij 



dridr2dr[dr'2 (l>*e^j^mj^ {r'lWejmj {r2)G{r[, r^, n, r2)(f)ejmj {ri)<Pe^j^mj^ {r2) 



_ ( i £h L] ( £ £h L' ] . I- I- . 

■Y^ 5^ jhSVlI'{-iy+^''--' \ 1/2 1/2 S W 1/2 1/2 S \ 

h JMLL'S [ j jh J ) [ j jh J ]■'■'■' ■' 
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x(-l) 
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In this case we define R\ and si as 



1 



(15) 



Introducing the same simplifying approximations as in 
the direct term, we obtain 



4(4- 



j j dkdk'ji 



X 



IQ2 
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(16) 



where Q2 and Q2 are defined by k and fe' by 



Q2 = - (^1 
Q2 = fl + 



2toi 

1712 ' 



m2 
2toi 



fe', 



2m 1 



fe' - ^fe. 

2mi 



(17) 
(18) 



These Q2 and Q'2 are obtained by changing the sign of 
fe' in Qi and Q[ of Eqs. (11) and (12). It is easy to sec 
that the difference of the expressions of Ie and Id is only 
the factor {-l)^+i. Thus, recovering the isospin degrees 
of freedom, we obtain 



+ = i(i^E/^^i \MRi)\'jH (1 + ^)' (^// dkdk'M\k + k'\R,)jlrfiM\)j^^MQ'. 



I) 



xP,,(cosQ^g2) ^ (1 + (-l)^+«+^-+V2-r)(/5MBl/2i^|TMB +i;,)2jgG^f ^(fc,fc')P,(cosfefe'). (19) 

qJgS 
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FIG. 1: Real and imaginary parts 
of single-particle potentials for 
and A in symmetric nuclear mat- 
ter at various Fermi momenta, 
0.75 fm-i <kF < 1.45 fm'^ 
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k [fm"'' 



In the above expression. Ib is the isospin of the baryon 
for which the s.p. potential is considered, and Mb is its z- 
component. The index denotes the proton or neutron 
in the target nucleus. 

The Eq. (19) defines a s.p. potential Ub{R) to give 

Id + Ie = 47r J R^dR \(j)e,{R)\^UBiR) (20) 

It is noted that the potential Ub{R) does not have £- and 
j-dependences due to approximations introduced in the 
derivation. 



III. SINGLE-PARTICLE POTENTIALS IN 
SYMMETRIC NUCLEAR MATTER 

Before discussing baryon s.p. potentials in finite nu- 
clei, we show s.p. potentials in nuclear matter at various 



Fermi momenta, 0.75 < kp ^ 1.45 fm"-'^, with the quark- 
model potential fss2 [l^. This potential is defined as a 
Born kernel of the RGM description of the interaction 
between the three-quark clusters. We use the energy- 
independent renormalized version of the fss2 potential 
[igj . The details of the G-matrix calculation for hyper- 
ons in nuclear matter are reported in ref. It has 

been known that the LOBT saturation curve in ordi- 
nary nuclear matter does not reproduce the empirical 
saturation property. Although the curve obtained by the 
potential fss2 with the continuous choice for intermedi- 
ate spectra almost goes through the empirical satura- 
tion point of kp = 1.35 fm"! and E/A = -16.5 MeV, 
the energy minimum appears at kp = 1.7 fm^^ and 
E/A = -20 MeV. Nevertheless the LOBT calculation 
provides a useful starting point and meaningful informa- 
tion for the baryon s.p. potentials in nuclear medium in 
the microscopic studies based on the bare baryon-baryon 
interactions. Missing effects in the LOBT , such as con- 
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FIG. 2: Real and imaginary parts 
of single-particle potentials for E 
and H in symmetric nuclear mat- 
ter at various Fermi momenta, 
0.75 fm-^ <kF < 1.45 fm'^ 



tributions from higher order diagrams and three-body 
forces are now semi-quantitatively understood in the nu- 
cleon sector [22| . 

Real and imaginary parts of the calculated s.p. po- 
tentials for A'^, A, S, and S in symmetric nuclear matter 
are shown in Figs. 1 and 2 as a function of the momen- 
tum k. These are the results after the self-consistency for 
the starting energy a; being reached. The Fermi momen- 
tum kp Yo chosen approximately in a step of one tenth 
of the normal nucleon density po- These densities are 
used as the discretized points of the density in the local 
density approximation for considering finite nuclei. Be- 
low kp = 0.75 fm^^, the nuclear matter G-matrix 
calculation becomes unstable due to the appearance of a 
bound state in the '^S\ channel. Because we expect little 
relevance of this phenomenon to ground states of finite 
nuclei, the instability is not inspected further. In the 
case that the effective interactions at low density below 



kp ~ 0.75 fm ^ is needed in the local density approxima- 
tion for finite nuclei, we use the G-matrices at kp = 0.75 
fm-i. 

As for the nucleon, the result is very similar to those by 
other realistic NN potentials. The depth of the s.p. po- 
tential at the normal density is considerably larger than 
the magnitude of the standard Woods-Saxon potential, 
which is 50 ^ 60 MeV. It has been known that the 
rearrangement potential reduces the strength by about 
10 ~ 20 MeV. 

We are concerned mainly with the prediction of the 
quark-model potential fss2 [l^ for hyperon s.p. poten- 
tials. The strength of the attractive A s.p. potential in 
normal nuclear matter shown in Fig. 1 is almost 45 MeV, 
which is again larger than the empirically known value 
of around 30 MeV. At the low density the potential is 
shallower. However, as will be shown in the next sec- 
tion, the depth of the calculated A s.p. potential in finite 
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nuclei, taking into account the finite geometry and the 
effects of the non-diagonal properties of the G, seems to 
be dictated by the potential depth at the normal nuclear 
density. We can expect that the rearrangement effects 
give rise to a repulsive contribution of the order of 10 
MeV for the A mainly through the energy dependence of 
the G-matrix. 

Nuclear matter calculations using the early version of 
the Kyoto-Niigata SUe quark- model potential, FSS, pre- 
dicted a repulsive S s.p. potential [20|- Results shown 
in Fig. 2 are obtained by the most recent quark-model 
parameterization, fss2 [1^. The E potential at fc = is 
definitely repulsive of about 15 MeV at normal density. 
This repulsion chiefly comes from the strong repulsive 
contribution of the '^S'l state in the isospin T = I chan- 
nel, which is naturally predicted by the quark-model as 
the consequence of the Pauli principle on the quark level. 
The interaction in the ^Sq with T = ^ channel is also 
repulsive. These repulsive contributions overwhelm the 
attractive contributions from the ^5*0 with T = | and the 
•^5*1 with T = i channels. The width T of the E hyperon 
in nuclear medium is related to the imaginary strength 
of the s.p. potential by T{k) = — 23[/(fc). r(0) is seen in 
Fig. 2 to be more than 30 MeV at normal density. 

The S s.p. potential in symmetric nuclear matter pre- 
dicted by fss2 is weakly attractive as is shown in Fig. 2. 
As the momentum increases, the magnitude of the at- 
traction is seen to become larger at the low momen- 
tum region of fc < 6 fm~^. The momentum dependence 
may be characterized by the effective mass. To obtain 
a rough estimation of it, we parameterize the potential 
by U^^^(k) ~ ak^ + b. In this case the effective mass at 
fc = is obtained by 




FIG. 3: Point nucleon density distributions pt{r) obtained by 
DDHF wave functions of the G-0 force 0|. 



2m~a 



(21) 




Calculated s.p. potentials give m~/ms ^ 1.1 at kp = 
1.35 fm-i and m^/ms ~ 1.05 at kp =1.07 fm^i. 

The S~p elastic and inelastic cross-section measure- 
ments at low energy by Ahn et al. [23l | indicate that 
the width of a S s.p. state in nuclear medium is F 3 
MeV. Although it is uncertain at which energy this num- 
ber should be compared with the calculated imaginary 
strength, the small imaginary strength of the S s.p. po- 
tential at the low momentum region given in Fig. 2 is in 
accord with the empirical small width of the S in nuclear 
medium. 

It is encouraging to see that the results for E and S 
hyperons agree at least qualitatively with empirical indi- 
cations so far obtained. 



IV. RESULTS IN FINITE NUCLEI 

We apply the calculational method presented in Sect. 
II to from light to medium-heavy nuclei: ^^C, ^^O, ^*Si, 
""^Ca, ^^Fe, and ^°Zr. Nucleon density distributions are 



FIG. 4: Real part of localized single-particle potentials for A'^, 
A, E, and H in ^^C with the quark-model potential fss2 
for the octet baryon-baryon interactions. 



prepared by density-dependent Hartree-Fock calculations 
using the Campi-Sprung G-0 force Q. Profiles of the 
point nucleon density distribution pt(r) which is a sum 
of the neutron and proton densities are shown in Fig. 3. 

Nuclear matter G-matrices are used in finite nuclei by 
the local density approximation. At the position Ri 
where the s.p. potential is evaluated the local Fermi 
momentum is defined by the correspondence kp[Ri) = 

, 2 , J 1/3 

^^pt{Ri) ■ The G-matrices calculated in nuclear 

matter with this Fermi momentum are used in Eq. (19). 
In actual calculations, G-matrix calculation is carried 
out only for the Fermi momenta shown in Figs. 1 and 
2. At each position the Fermi momentum which 
is closest to kp{Ri) among these twelve values is cho- 
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FIG. 5: Same as in Fig. 4, but for ^^O. 



FIG. 7: Same as in Fig. 4, but for '^''Ca. 





FIG. 6: Same as in Fig. 4, but for ^^Si. 



FIG. 8: Same as in Fig. 4, but for ^'^Fe. 



sen. As explained in Sect. II, for small densities below 
kp — 0.75 fm^^, namely the total density pt = 0.028 
fm^'^, we always use kp = 0.75 fm~^. In homogeneous 
matter the s.p. potential is determined by the matrix 
element with the zero momentum transfer, namely diag- 
onal {k' = k) components of G{k' ,k; K,uj). In finite 
nuclei non-diagonal components of the G matrices also 
contribute to the s.p. potential. 

The starting-energy dependence of the G-matrix plays 
an important role in the LOBT. The prescription of the 
starting-energy as the sum of s.p. energies of the two 
baryons considered means that certain higher-order dia- 
grams are included. Hence the self-consistency between 
the s.p. energy which is defined by the G-matrix and 
the G-matrix which depends on the starting-energy is re- 



quired. Calculations in nuclear matter shown in Sect. 
II are the results with this consistency achieved. In 
the case that the G-matrix in nuclear matter is applied 
to a finite nucleus, however, there is no simple way to 
treat the starting-energy dependence. We introduce an 
ad hoc prescription to use an s.p. potential value at 
the median momentum 2^^^^kp: uj = 2niM H — —k^ + 

+ UN{2-^'^kp) + UNi^-^'^kp) for the nucleon 



and oj = TOy 



niN 



7\ ft- 



2(mi-4-*TiAf) 



C/y(0) -I- UN{2~^'^kp) for the hyperons. 

The results with the SUq quark-model potential fss2 
in the energy-renormalized form are shown in 
Figs. 4~9. The charge state of the baryon specified by 
Mb in Eq. (19) is set to be Mb = —Ib- Comments on 
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8 



r [fm] 



FIG. 9: Same as in Fig. 4, but for ■'"Zr. 



the calculated s.p. potential of each baryon are given in 
the following. 



A. Nucleon s.p. potential 

The shape of the calculated nucleon potential is seen to 
follow the density distribution, and the depth is 80 ~ 90 
MeV which corresponds to the s.p. potential in nu- 
clear matter at the normal density. It is well known 
that the straightforward application of the LOBT start- 
ing from realistic NN interactions overestimate the at- 
tractive nucleon-nuclcus potential. To compare the cal- 
culated potential with the empirical one, we need to in- 
clude the so-called rearrangement potential. The repul- 
sive strength is known to be 10 ~ 20 MeV Q. If this 
contribution is taken into account, the resulting poten- 
tial becomes closer to the phenomenological potential of 
the Woods-Saxon form. 

It is noted that it is still a remaining problem in nu- 
clear physics to understand nuclear bulk properties in a 
fully microscopic way on the basis of the realistic inter- 
actions including higher-order correlations, three-body 
forces, and other possible medium effects. 



B. A hyperon s.p. potential 

The AA^ ^5*0 state has a similar character to the NN 
^So state in the spin- flavor SUg symmetry, although there 
is small admixture of a completely quark Pauli forbid- 
den component. Similarly, the interaction in the AA'^ ^Si 
channel resembles that in the NN channel with a 
smaller magnitude by a factor of l/V^, although there 
is an important difference that the pion exchange is ab- 
sent. Thus it is probable that the A-nucleus potential 



is attractive, about a half of the A''-nucleus potential in 
magnitude. Looking at the density-dependence of the 
A s.p. potential in nuclear matter, we should expect a 
similar rearrangement potential as in the nucleon case. 
The addition of the hyperon to nuclear medium does not 
change directly the nucleon density and hence the nu- 
cleonic Pauli effect. The rearrangement effect for the 
A hyperon originates, in the LOBT, mainly from the 
energy-dependence of the NN and AN G-matrices. If 
we assume a repulsive rearrangement potential of the or- 
der of 10 MeV, calculated results shown in Figs. 4^9 
correspond well to the empirical A-nucleus potential in 
the Woods-Saxon form with the depth of about 30 MeV. 



C. E hyperon s.p. potential 

The experimental information has been limited for the 
S s.p. potential in nuclei. Because the S state in a nu- 
cleus is expected to have a large width due to the strong 
UN AN conversion process, it is unlikely to observe 
clear peak structure in the E-formation spectra. Never- 
theless, results from the early experiments of {K, vr) in- 
clusive spectra measured at CERN were interpreted 
as indicating that the E-nucleus potential is moderately 
attractive. The discovery of the bound ■^He-I-I]*' and 
'^II-|-S+ systems 1251 and the theoretical consideration 
by Harada et al. [2y| showed that the attraction in the 
UN T — 1/2 channel should be attractive enough. An- 
other important source of the S-nucleus interaction is the 
energy shift and the width of atomic orbits extracted 
from the X-ray data. Batty, Friedman, and Gal [23| an- 
alyzed the data to conclude that the E-nucleus potential 
changes its sign toward higher density region in a nu- 
cleus from the attractive potential at the surface region. 
Dabrowski [28| aiialyzed the BNL experiment of {K, tt) 
spectrum on ^Be [23| in a plane wave model and conjec- 
tured that the S potential is repulsive of the order of 20 
MeV. Recent experimental data with better accuracy of 
(tt^, X+) inclusive spectra measured at KEK [s^l was re- 
ported to suggest that the E-nucleus potential is strongly 
repulsive, the strength being more than 100 MeV. Sev- 
eral theoretical analyses carried out later [3l|, [32| con- 
firmed the repulsive nature of the E s.p. potential, but 
the height may be a few 10 MeV. 

On the theoretical side, EA^ interaction models admit 
large uncertainties. Most parameter sets of the Nijmegen 
hyperon-nucleon OBEP potential [1, H, [l^, [13 predict an 
attractive E-nucleus potential in nuclear matter, except 
the model F [sl, [H, [1^1 . On the contrary, the strong 
repulsive character in the T = 3/2 ■^5jchannel is inher- 
ent in the quark-model description [3, [s^. Thus the 
Kyoto-Niigata quark-model potential predicts an overall 
repulsive E s.p. potential in nuclear matter. Results 
given in Figs. 4~9 show the consequence of this property 
to finite nuclei. At higher density region inside a nucleus 
the E-nucleus potential is repulsive of 10~20 MeV. The 
overall repulsive nature of the E-nucleus potential has 
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been deduced by the analyses of {K~ ,Tr~^) S~ formation 
inclusive spectra [s^, HH, HI] . Beyond the surface region 
the potential becomes attractive. It is interesting to see 
that the radial dependence indicated by the analyses of 
E atomic data 27] is actually derived by the microscopic 
calculation using the quark-model bare interaction with 
no phenomenological adjustment. 

It is remarked that in the present evaluations we ap- 
ply G-matrices in symmetric nuclear matter to finite nu- 
clei without separating neutron and proton density dis- 
tributions. However, in heavy nuclei, e.g. ^"Zr in our 
calculations, in which neutron and proton distributions 
are visibly different, we should take care of the isospin- 
dependence when solving the Bethe-Goldstone equation. 
In that case, the repulsive contribution to the E" s.p. 
potential from the T = 3/2 ^Si channel becomes more 
predominant. 



D. H hyperon s.p. potential 

For the 2 s.p. potential, the experimental informa- 
tion has been more scarce than the E. The S hyperon is 
formed in (K~,K~^) reaction on nuclei with small pro- 
duction rates. There has been no concrete evidence of the 
S hypernuclear bound state. The existing experimental 
data of the S formation spectrum 38, 39] has sug- 
gested that the S feels attractive potential in nuclei, the 
depth of which is not so large, 10 ~ 20 MeV. Our results 
shown in Figs. 4~9 with the quark-model potential fss2 
[m show that the S-nucleus potential is weakly attrac- 
tive at and beyond the nuclear surface region, which is 
similar to the E-nucleus potential. Toward the inside of 
the nucleus the H s.p. potential tends to be repulsive 
and oscillates around zero with the magnitude of about 
10 MeV. It is not possible to simulate the potential shape 
by a single Woods-Saxon form. No S hypernuclear bound 
state is expected from such a weak potential. The situa- 
tion does not change even if the actual potential strength 
differs by a factor of 2 or so. In that case the level shift 
of the S~ atomic orbit should be a valuable source of 
the information about the S-nucleus interaction. This 
subject is addressed in the next section. 

The evaluated SiV G-matrices include full baryon- 
channel couplings, namely the possible S7V-AA-EE or 
SA^-AE-EE couplings. It is helpful to use equivalent in- 
teractions in low-momentum space to check the charac- 
ter of the EN interaction and the effect of the baryon- 
channel coupling in each spin and isospin state. Inspect- 
ing the matrix elements in ref. 40], we see that the EN 
effective interaction from the fss2 in the T ~ 1 channels 
both in ^5*0 and ^Si are repulsive. The T = ^5*1 in- 
teraction is very weak, and the T = ^5*0 interaction is 
attractive for which the S7V-AA-EE coupling is respon- 
sible. It turns out that the net s-wave contribution is 
small and thus the attractive p-wave contribution plays 
an important role to make the E s.p. potential to be 
attractive at the surface region. 



TABLE I: The strength and geometry parameters of the 
Woods-Saxon form fi(r) = U° /[I + exp((r — ro,i)/«0] fit- 
ted to the real part as well as the imaginary part of E and H 
s.p. potentials UB{r) in ^*Si and ^^Fe. 
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-25.94 4.179 
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-f57.43 3.049 


0.7860 
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0.2948 
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-3.746 6.035 


0.5655 


-65.78 


4.955 


0.8902 
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22.19 4.031 


0.3990 


+37.65 


5.486 


0.8867 




3 


32.22 1.553 0.5527 +16.77 


-0.4167 


1.163 




1 


-2.232 6.367 


0.2389 


-15.92 


5.845 


0.6820 




2 


+8.295 3.796 0.2597 +14.64 


5.947 


0.6728 




3 


+20.67 1.674 


0.4887 


+1.290 


0.7581 


0.1709 



V. ENERGY SHIFT AND WIDTH OF ATOMIC 
ORBIT 

The level shift and the width of E^ atomic orbits are 
a valuable source of the information on the E-nucleus 
strong interaction. The analyses by Batty, Friedman, and 
Gal [23 indicated that the E-nucleus potential is attrac- 
tive at the surface region, but at higher density region in 
a nucleus the potential turns to be repulsive. The radial 
dependence of the calculated E~ s.p. potential shown 
in the previous section agrees with this. Therefore it is 
instructive to explicitly evaluate the energy shift and the 
width of E~ atomic orbits with the calculated potential. 
As will be shown below, the result is consistent with the 
experimental data. This indicates that the microscopic 
calculation with the quark-model fss2 is reliable in the 
EiV channel. Thus, it is interesting to extend the level 
shift calculation to E~ atomic orbits. The experimen- 
tal data should be available in near future, because the 
first measurement of E~ atomic X rays from Fe target 
is proposed [13] to be performed at J- PARC. The theo- 
retical prediction provides a guiding information for this 
experiment. 

We consider ^®Si and ^^Fe for explicit evaluations of 
the level shift of the atomic orbit. We first fit the shape 
of the calculated s.p. potential using the Woods-Saxon 
form. For ^^Si a sum of three Woods-Saxon shapes is used 
and for ^^Fe a sum of two Woods-Saxon shapes and one 
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derivative of the Woods-Saxon shape is assigned. Param- 
eters are given in Table I. It is noted that the imaginary 
parts are also given to illustrate the order of the magni- 
tude of the absorptive strength, intending to demonstrate 
that the S imaginary potential is about one order of mag- 
nitude smaller than the E one. However, actual num- 
bers should not be taken very seriously because nuclear 
matter calculations tend to overestimate the imaginary 
strength as the calculations of nucleon optical model 
potential indicates. In addition, the prescription to use 
the G-matrices at kp = 0.75 fm^^ for all the densities 
below kp = 0.75 fm^^ probably leads to the overestima- 
tion of the imaginary strength at the surface region. It is 
also remarked that localized imaginary potential through 
the zero-momentum Wigner transformation may become 
positive at some points. 



A. E" 

Results of the level shift AE — E ~ Ec and the width 
r = —2QE for the S~ /- and (7-atomic levels on ^^Si and 
the T," g- and /i-atomic levels on ^^Fe are given in Table 
II, where Ec stands for the Coulomb bound state energy 
without the I]~-nucleus strong interaction. When the 
real part of the E s.p. potential is taken into account, 
the energy of the n = 4,1 = 3 orbit on ^*Si is shifted 
downward by 222 eV. To investigate the contribution of 
the absorptive effect, we do not use the calculated poten- 
tial given in Table I. The imaginary potential is rather 
strong as explained above. The magnitude of the level 
shift and the width depends non-linearly on the strength 
of the imaginary potential. Hence we use an phenomeno- 
logical imaginary potential to discuss the level shift of the 
atomic orbits of the E~. We add an imaginary potential 
in the single Woods-Saxon form used in ref . [27| , namely 
To = 1.1 X 28^/3 and a = 0.67 with the depth of -9 
MeV. In that case, we obtain AE = 208 eV and T = 249 
eV, which well correspond to the experimental values of 
AEe:,p = 159 ± 36 eV and Te^p = 220 ±110 eV. This 
result indicates that the real part of the E s.p. potential 
calculated microscopically in the LOST starting from the 
two-body quark-model potential fss2 is reasonable. 



B. S" 

Observing that our calculated E-nucleus potential 
gives a reliable result for the shift of the E~ atomic level, 
it is interesting to proceed to S~ atoms without any ad- 
justment. In ref. [43] Batty, Friedman, and Gal esti- 
mated the level shift and the width of S~ atoms for the 
S^-nucleus potential having an attraction of the depth of 
15-20 MeV with an imaginary strength of 1-3 MeV in the 
shape roughly following a nuclear density distribution. 
Although the S-nucleus potential is weakly repulsive in- 
side the nucleus in our calculation, the attractive strength 
at the surface region is found to be comparable to that 



TABLE IL The energy shift AE ^ E - Ec and the width 
r = -2QE of the E" atomic orbits in ^**Si and ^''Pe, using 
the parameterized E s.p. potential given in Table I for the 
real part. The imaginary potential is given in a Woods-Saxon 
form with the strength Wo = —9 MeV and the geometry 
parameters ro = 1.1^4^'''^ fm and a — 0.67 fm. Entry numbers 
are in eV. 
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-943 
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real part only 




-11 




real -|- imaginary 




-11 


8.3 



of the Woods-Saxon potential with the depth of 10 ^ 20 
MeV and the geometry parameters of ro = I.IA^/'^ fm 
and a=0.67 fm. Results for the /- and g-orbits in ^*Si 
and the g- and /i-orbits in ^^Fe are given in Table HI, 
together with those of the Woods-Saxon potential as a 
reference. Because the atomic level shift is insensitive 
to the short range part of the strong interaction poten- 
tial, our potential parameterized as shown in Table I, 
predicts a similar magnitude of the level shift of the ref- 
erence Woods-Saxon potential. Because the S imaginary 
potential obtained by the G-matrix is small, we directly 
use it for the estimation of the width F, though it is likely 
to overestimate the absorptive effect as in the case of E~ . 
As is seen in Table HI, the width is of the order of a few 
hundred eV for the g-orbit in ^^Fe and the energy shift 
is hardly affected. It is needed to see what order of the 
magnitude is detected for the S~ n = 5,£ = A level in 
^^Fe in the future experiment prepared at J-PARC (43| . 

VI. SUMMARY 

In order to examine the prediction of single-particle 
properties of all the octet baryons in nuclear medium, 
especially S hyperon, by the recently developed quark- 
model baryon-baryon interactions, we have evaluated lo- 
calized s.p. potentials in finite nuclei by folding G- 
matrices in nuclear matter with respect to nucleon s.p. 
wave functions in the scheme of the local density approx- 
imation. Introducing a spin-average approximation and 
a zero-momentum Wigner transformation, the non-local 
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TABLE III: The energy shift AE ^ E - Ec and the width 
r — —2'SE of the H~ atomic orbits in ^**Si and ^''Fe, using 
the parameterized H s.p. potential given in Table I both for 
the real and imaginary parts. As a reference, results obtained 
by the complex Woods-Saxon potential with the strength of 
Uo = -14 - 3i MeV, ro = I.IA^^^ fm, and a=0.67 fm are 
included. Entry numbers are in eV. 
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baryon s.p. potential calculated in momentum space is 
reduced to a local potential in coordinate space. The fi- 
nal expression is feasible for numerical calculations, in the 
case that the nucleon density distribution is discretized. 
Adopting about one tenth of the normal nuclear den- 
sity as the interval of the discretization, we have car- 
ried out calculations in ^^c, ^^O, ^^Si, '^°Ca, ^^Fe, and 
^"Zr for each octet baryon; N, A, E, and S. This is the 
first comprehensive evaluations of the s.p. potentials of 
all the octet baryons in finite nuclei, starting from the 
baryon-baryon bare interactions. These microscopic cal- 
culations of octet baryon s.p. potentials in finite nuclei 
are meaningful to elucidate the character of the theo- 
retical model of the octet baryon-baryon interactions by 
comparing them with empirical s.p. potentials, which are 
not yet available for the E and S hyperons. 

We use the most recent quark-model potential, fss2 
[l8| as bare baryon-baryon interactions. The energy- 
dependence in the original form of the quark-model po- 
tential is eliminated by the renormalization procedure. 
The NN sector of this potential describes scattering data 
as accurately as other modern realistic interaction. Cal- 
culated nucleon s.p. potentials in nuclear medium are 
found to be similar to those obtained in the LOBT frame- 
work with other potentials. The A7V interaction is under 
control to a certain extent by the experimental data of 
A hypernuclei. The fss2 gives similar A s.p. potentials 
in nuclear medium to those of the Nijmegen OBEP po- 



tential. It was shown [i^ in fact that the fss2 and the 
Nijmegen NSC97f actually have very similar matrix el- 
ements of the equivalent interaction in low-momentum 
space. 

The extension to the 'SN channel and further to the 
SN sector of the strangeness S* = — 2 has many ambigui- 
ties because of scarcity of experimental information. It is 
necessary to rely on the theoretical framework as reliable 
as possible to construct these baryon-baryon interactions. 
At present, the SUg quark-model fss2 is more predictive 
than the OBEP model, in the sense that the potential 
parameters are uniquely given in contrast to various sets 
of parameters presented by the Nijmegen group. Thus 
we focus our attention in this paper on the fss2 as the 
input TiN and two-body interactions. The compari- 
son with the results by other potential models including 
the parameterization based on the effective chiral field 
approach [H, is an interesting future subject. 

Properties of the calculated S s.p. potential in finite 
nuclei are found to agree well with the experimental ev- 
idence so far obtained. The potential is repulsive of the 
order of 10-20 MeV, which is necessary to accounts for 
the (7r^,A'+) E^ formation spectra on nuclei [3l|, [33 |. 
while it should be attractive at the surface region, as the 
atomic level shifts indicate [131 ■ The attraction obtained 
by the present microscopic calculation can reproduce the 
empirical energy shifts of the S~ atomic levels. To de- 
termine more precisely the shape and the strength of the 
E-nucleus potential including its isospin-dependence and 
study the relation to the underlying bare interaction, we 
need further experimental data. 

On the basis of these observations, it is interesting to 
consider the EN sector. The imaginary part of the S 
s.p. potential is small in the quark model description, 
which is consistent with the empirical estimation based 
on the S~p scattering cross-section at low energy p3| . 
Therefore, if the real part is sufficiently attractive to sus- 
tain hypernuclear bound states, we can expect clear spike 
structure in the (A'^, K^) inclusive spectra on nuclei. At 
the surface region the potential has similar attraction to 
the E potential. However, the potential does not have 
a familiar shape simulated well by a Woods-Saxon form. 
Inside the nucleus the potential fluctuates around zero, 
reflecting the fluctuation of the nucleon density distribu- 
tion. The quark model potential fss2 implies that the net 
s-wave contribution to the S s.p. potential is small and 
the net weakly attractive contribution from the p- waves 
is relatively important. Such a potential does not sup- 
port S nuclear bound states. In that case, measurements 
of the energy shift and width of the S atomic level be- 
come an invaluable source of the information about the 
S-nucleus strong interaction. Our S s.p. potential calcu- 
lated by the fss2 suggests that the negative energy shift 
is the order of a few hundred eV for the 4/-orbit in ^^Si 
and the order of 1 keV for Sg-orbit in ^^Fe. The existence 
or non-existence of S hypernuclear bound states and the 
energy shifts of the S~ atomic states will be clarified in 
near future by the experiments prepared in the J- PARC 
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project [47|, which should advance our understanding of 
the interactions in the S = —2 sector. Our calculational 
framework provide a useful method to link the baryon 
s.p. properties in nuclei with the two-body interactions 
between octet baryons. 
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